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Abstract. Recent HST determinations of the expansion's 
rate of the Universe (the Hubble constant, Hq) assumed 
that the Cepheid Period-Luminosity relation at V and I 
are independent of metallicity (Freedman, et al., 1996, 
Saha et al., 1996, Tanvir et al., 1995). The three groups 
obtain different vales for Hq. We note that most of this dis- 
crepancy stems from the asumption (by both groups) that 
the Period-Luminosity relation is independent of metallic- 
ity. We come to this conclusion as a result of our study 
of the Period-Luminosity relation of 481 Cepheids with 3 
millions two colour measurements in the Large Magellanic 
Cloud and the Small Magellanic Cloud obtained as a by- 
product of the EROS microlensing survey. We find that 
the derived interstellar absorption corrections are partic- 
ularly sensitive to the metallicity and when our result is 
applied to recent estimates based on HST Cepheids obser- 
vations it makes the Iow-TJq values higher and the high-_ffo 
value lower, bringing those discrepant estimates into agre- 
ment around Hq « 70 kms~^Mpc~^. 
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1. Introduction 

The value of the Hubble constant Hq fixes the fundamen- 
tal physical scale and time scale of the Universe, there- 
fore by inference its age. Dichotomous results are obtained 
in discordant ranges, even allowing for their internal er- 
ror estimates. The results are dependent upon the choice 
of the cosmological distance ladder used (Fukugita et al. , 
1993, Freedman et al., 1996, Saha et al., 1996). However, 
both the longer distance scale {Hq « 50 kms~^Mpc~^), 
and the shorter distance scale {Hq « 80 kms~^Mpc~^) 
rely strongly on the Cepheid period-luminosity relation 
calibrated in the Large Magellanic Cloud (LMC), and 
assumed to be independent of metallicity. The Cepheid 
Period-Luminosity (PL) relation is widely considered to 
be the most accurate way available to measure extragalac- 
tic distances on small scale. The Cepheids are relatively 
young, bright periodic variable stars which pulsational pe- 
riods are strongly correlated with their luminosity. Their 
distance is therefore determined via a comparison to a 
calibrated PL relation. Cepheids are observed in distant 
galaxies, up to the Virgo cluster, in order to calibrate 
secondary distance indicators that connect to the Hubble 
flow. 

The theory (Stothers, 1988, Stiff 1990, Chiosi et al., 
1993, Stiff 1995) predicts a small abundance ( helium frac- 
tion and heavy elements) effect on the PL zero point due 
to : (i) theory of stellar pulsation through the dependence 
of period on mass and radius, (ii) theory of stellar evolu- 
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tion through the mass luminosity relation, (iii) theory of 
stellar atmosphere through the relation between effective 
temperature, absolute magnitude in bandpasses and bolo- 
metric correction. Disentangling the effects of metallicity 
and reddening has a venerable history (Feast 1991). Vari- 
ous observational studies since the 70s, have looked for and 
found a color shift between LMC and SMC Cepheids (e.g. 
Gascoigne 1974, Martin et al. 1979). Caldwell & Coul- 
son (1986) and more recently Laney & Stobie (1994) used 
theoretical and empirical relations, as well as individual 
reddenings from color-color diagrams to derive PL rela- 
tions adjusted for abundance differences. Their approach 
is different from ours (and the technique used by the HST 
teams); a full comparison will be model dependent (flux 
redistribution models) and we have not attempted to do 
that here. A limited comparison to the VJHK Table 6 of 
Laney & Stobie (1994) shows that we confirm the sense of 
their results. 

An empirical test (Preedman and Madore, 1990, 

Gould, 1994) in three fields of M31 with 36 Cepheids and 
152 BVRI measurements have led to ambiguous results. 
As part of the HST key project on extragalactic distance 
scale observations, a further check is planned by obser- 
vations of two fields in MlOl. Currently, all HST studies 
assume that the Cepheid PL relation has no metallicity de- 
pendence (Freedman et al., 1994a, Freedman et al., 1994b, 
Sandage et al., 1994, Tanvir et al., 1995, Kennicutt et al., 
1995). 

2. EROS observations. 

We used a new data set of 481 Cepheids and 3 millions 
measurements in two colours obtained in the LMC and the 
SMC to derive the dependence of the optical period lumi- 
nosity relations on metallicity. Observations were obtained 
as a by-product of the EROS (Experience de Recherche 
d'Objets Sombres, Aubourg et al., 1993) microlensing sur- 
vey. 

Observations have been obtained from ESO La Silla 
using a 0.4m f/lO reflecting telescope and a 2 x 8 mo- 
saic of CCDs (Arnaud et al., 1994ab). During the ^ 100 
nights of 1991-1992 campaign (Beaulieu et al., 1995, Gri- 
son et al., 1995), about 2000 images have been obtained in 
the bar of the LMC in a blue and a red bandpass {Be and 
Re)i whereas about 6000 images have been obtained in a 
crowded field of the SMC with a near pair of very similar 
fihers {Be2 and Re2) during the - 200 nights of 1994- 
1995 season (Beaulieu et al., 1996). An accurate and re- 
liable photometric transformation have been constructed 
between the two systems. Hence we present a differen- 
tial analysis that relies on no external zero-point and thus 
completed within the {Be and Re) system. Since we apply 
our result to HST V , I observations, we make use of the 
fact that the net transmission of the Be band after convo- 
lution with the CCD reponse is much closer to Johnson V 
than to Johnson B. The transformation is similar to that 



between the broader HST F555W filter and V. The Re 

band is between Cousin R and /. The Be — Re color trans- 
forms wen to V-I : V-I = 1.02{Be-Re), o = 0.02mag. 

The two data sets were searched for Cepheids with 
Fourier analysis (Grison 1994) and phase dispersion min- 
imisation techniques (Schwarzenberg-Czerny 1989). Stars 
affected by blending effects have been excluded of the sam- 
ple. The high quality, excellent phase coverage light curves 
were Fourier analysed to separate between the classical 
Cepheids that pulsate in the fundamental mode, and the 
so called s-Cepheids that pulsate in the first overtone, and 
therefore follow a different period luminosity relation. In 
the LMC we keep 51 fundamental pulsators and 27 first 
overtone pulsators, and 264 fundamental pulsators and 
141 first overtone pulsators in the SMC. So, we have two 
complete samples of Cepheids that fill densely the period- 
magnitude-color [PLC] space, with known difference in 
metallicity IS.\Fel H\i,mc-Smc = 0.35 (Spite and Spite, 
1991, Luck and Lambert 1992). 

3. The method 

We compare the two samples in the PLC manifold to de- 
rive 2 independant sources of difference — distance and 
extinction, and to search for a third one — metallicity. 
The method will be applied independently to both funda- 
mental and first overtone pulsators, and is a minimiza- 
tion fit to a model. Details of the method are presented 
elsewhere (Sasselov et al., 1996). The model is built on 
the basis of the technique used for determining Cepheid 
distances to galaxies by the HST Key Project (Preed- 
man and Madore, 1990). It is as follows: the LMC PL 
relation is used to slide the observed PL relation against 
it; the Galactic extinction law is applied in deriving the 
reddening (implicitly assuming no difference between the 
ensemble colors of the Cepheids in LMC and the observed 
galaxy); a true distance modulus is derived by adopting a 
LMC distance and correcting for extinction (using the red- 
dening derived above with adopted LMC reddening E(B- 
V)=0.10). Other assumptions in the technique (and our 
model) are: constant PL slope with metallicity (confirmed 
by our observations) and no depth dispersion in the LMC 
Cepheid sample (which is only from the LMC bar). In the 
application we use as constrains the independently derived 
foreground reddening of LMC (0.06) and SMC (0.05), and 
the line of nodes for SMC from Caldwell & Coulson (1986) 
to derive the EROS SMC sample depth dispersion. We 
fit for different fixed values of the extinction parameter 
Rv- The model describes the case considered by Stothers 
(1988) in his equation (26). 

We model the Be, Re data in the PLC space tak- 
ing into account the high degree of correlation between 
the measurements following the formalism introduced by 
Gould (1994). Unlike Gould, we solve for a wavelength 
dependance of a metallicity effect. Our model has 12 pa- 
rameters, which are : linear fits of the corelations in the 
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PLC projections (slopes /?i,6i and the zero points ai,ai), 
the distance difference 71, the relative reddening differ- 
ence 72, and the metallicity terms 7I {i=l,2 are the two 
passbands, Be and Re)- The model parameters are esti- 
mated by minimizing the mean magnitude residuals, Xj.p, 
where Qi^p is the observed mean magnitude in the i band 
for a Cepheid with period Pp. For example, the residuals 
for the SMC sample are given as: 

^Ip = Qi,P - (a< + PilogPp + 71 + j2Ri + 73). (1) 

^tp = Qt,p-[ai + hiQi,p-Q2,p + j2iR2-Ri) 

+73 - 73') + 71 + l2Ri + lU , (2) 

with p = 1, N{n) being the number of Ccpheids. The 
form of the covariance matrices of the data, the mini- 
mization and the iteration procedure are described in Sas- 
selov et al. (1996). 

4. The results 

We applied the global fitting procedure to the data-set of 
fundamental mode and first overtone mode pulsators in- 
dependently, and we obtained the same results. 

The two set of PL relation give, before correction for 
metallicity and reddening, the same distance modulus dif- 
ference of 6{LMC - SMC) = 0.89 ± 0.05mag. The other 
model parameters have the following values for the fun- 
damental pulsators : ai,(3i = 17.61 ± 0.035, -2.72 ± 0.07; 
17.74 ±0.029, -2.95 ± 0.06; ai,6i = 14.49 ±0.13, -11.78 ± 
0.94; 14.36 ± 0.12 - 12.77 ± 0.90. 71 = 0.62 ± 0.04, 72 = 
0.01 ± 0.01, 7,^ = 0.06 ± 0.01, 7f = -0.01 ± 0.01. We 
adopted the same extinction law with Ry = 3.3 for the 
SMC. The correction due to the metallicity dependence of 
the inferred distance modulus by the described technique 
is then Sjj, = A/itrue - ^fJ-infered = -0.139tQ;ogmag. 

The sources of error are the photometric transforma- 
tion, between the two data sets (0.06 mag), the scatter 
in the PL relation (0.04 mag) and the uncertainties in 
the reddening and depth determination in the SMC (0.03 
mag). 

An alternative to the metallicity effect could be un- 
usual SMC extinction with a value of Ry > 5 and 
E{B — V) < 0.04. We consider this alternative to be 
very unlikely, because it requires that the extinction be 
less than the foreground to SMC and an unusually high 
Rv value unjustified by any other evidence. On the other 
hand, with Ry = 3.3, our mean reddening value is in good 
agreement with a number of other independent estimates 
(Bessell 1991). 

We derive the following metallicity dependence of red- 
dening corrected distance moduli of Cepheids: 

<5/z = (0.4410:1) log-^, 

^LMC 



where Z is the abundance of heavy elements (by mass) in 
the studied Cepheids and Zlmc=0.0085. The metallicity 
dependence is valid in the spectral region covered by the 
EROS filters, but can be applied to HST (V,I) work for the 
reasons presented in Sect. 2 above. The form of the metal- 
licity dependence used here is based on theoretical consid- 
erations and is a good assumption in the suggested range 
of appHcation of 0.001 <Z< 0.02 (Chiosi et al. 1993). 

Despite its relative weakness, this metallicity effect is 
significant to extragalactic distance measurements. This 
is especially true when the color shift due to metallicity 
73 — 7I is interpretated as reddening in the determination 
of the true distance modulus of a target galaxy (Preedman 
et al., 1994a, Freedman et al., 1994b, Sandage et al., 1994, 
Tanvir et al., 1995, Kennicutt et al., 1995). 



5. Effect on Hq determination 

Recent efforts to determine iJg from Cepheid distances 
have focused on: (1) the Virgo cluster galaxies (Freedman 
et al., 1994a), (2) the Leo I group (Tanvir et al., 1995), 
containing elliptical galaxies; and (3) the parent galaxies 
(Sandage et al., 1994) of supernovae of type la. All these 
studies use the same modern technique with the LMC as a 
base and all the same initial assumptions; we share them 
in our analysis. Yet they result in three different values of 
Hq, with hardly overlaping error bars. The abundances of 
metals in the young populations of the host galaxies for 
these studies differ and, apparently, in a systematic way. 
Therefore we propose that the metallicity dependence may 
be responsible for most of this discrepancy. 

With the first approach and HST VI photome- 
try of Cepheids in MlOO, they derive ilo = 80 ± 17 
km.s~i Mpc~^. For the metallicity of the Cepheids we 
adopt [Fe/H]=+0.1 (Z=0.02), derived from the abun- 
dances of 

H II regions (Zaritsky et al., 1994) and assuming a so- 
lar [0/Fe] ratio. With the uncertainty range of our metal- 
licity effect, this leads to Ho = 76 - 70 km.s"^ Mpc~^. 

With the second approach and HST VI photometry of 
Cepheids in M96 (Tanvir et al., 1995), they derive Hq = 
69 ± 8 km.s-i Mpc'^. Oey & Kennicutt (1991) give an 
abundance estimate for M96 of [Fe/H]=-0.02. This implies 
a small change in the Cepheid distance to M96 and the 
value of Ho, i?o = 64 - 66 km.s'^ Mpc''^. 

With the third approach and HST VI photometry of 
Cepheids in IC4182 and NGC5253 (Sandage et al., 1994), 
they derive i?o = 55 ± 8 km.s"^ Mpc~^. Abundances in 
NGC5253 have been measured and discussed (Pagel et al., 
1992); discussion of abundances in IC4182 has been also 
given by Saha et al., (1996) — we adopt [Fe/H]=-1.3 and 
caution on the large uncertainties. With this metallicity 
Hq = 59 — 68 km.s~i Mpc~^. Other systematic effects 
in the SN la method (Riess et al., 1995), bring the value 
of Hq even further up; to a value also favored by new 
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theoretical SN models (Hoflich and Khokhlov 1996). The 

new value of Hq = 58 ± 4 km.s"-'^ Mpc~^ (Sandagc ct 
al., 1996) includes two SNe la in spiral galaxies, which 
Cepheid metallicity is most likely similar to that of the 
LMC. Note, however, that the two SN (1981B and 1990N) 
are by 0.2—0.4 mags dimmer than the three SNs in IC4182 
and NGC5253. Therefore, the result we obtained above 
by accounting for the effect of metallicity on the Cepheid 
distances for each galaxy, remains virtually unchanged. 

The metallicity dependence we found from the 
LMC/SMC analysis brings all the derivations of Hq to 
good agreement. It is summarized in Fig. 2. 



LMC/SMC 



Ho (km/s/Mpc) 
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0.2 0.4 



Fig. 1. The colour magnitude diagram for the fundamental 
mode Cepheid pulsators in the SMC (diamond) and LMC 

(filled diamond) in the Be — Re colour system. The differ- 
ence vector represents the reddening line in this plane. The 
color difference between the two samples of Cepheids is due 
to their different metallicity. The SMC Cepheids are bluer and 
fainter at a given pulsation period. 
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